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SECTION"  I 


INTRODUCTION 

One  VLF  magnetic  loop  receiving  antenna  on  the  EC- 135  aircraft 
consists  of  two  orthogonal  conductor  coils  wound  around  a rectangular 
slab  made  up  of  two  arrays  of  thin  ferrite  rods.  The  rod  assembly  is 
mounted  flat  on  the  metallic  skin  of  the  fuselage.  Essentially  the  func- 
tion of  the  ferrite  is  to  concentrate  the  magnetic  flux  of  the  incoming 
VLF  signal.  The  receiving  antenna  is  completely  characterized  by 
(a)  the  input  im.pedance  Z.^  which  is  mainly  the  coil  self-inductance, 
and  (b)  the  open- circuit  voltage  which  is  proportional  to  the  maa- 

netic  flu.x.  These  two  quantities  can  be  calculated  by  solving  two  quasi- 
static boundary  value  problems  of  a rectangular  ferrite  slab  in  applied 
magnetic  fields. 

In  t.his  report  the  general  problem  of  a perm.eable  body  in  an  exter- 
nal quasistatic  magnetic  field  is  first  formulated  as  a Fredholm  integral 
equation  of  the  second  kind.  For  a uniform  external  field  this  integral 
equation  is  solved  approximately  by  the  method  of  averaging  functional 
corrections  [1].  The  fl'ux  enhancemient  factor  of  a highly-permeable 
rectangular  slab  is  determined  for  a wide  range  of  the  thic.kness-to- 
length  ratio.  The  same  approximation  method  is  also  used  to  calculate 
the  self- inductance  of  a thin- wire  current  loop  wound  around  ^he  mid- 


section  of  the  rectangular  slab. 


SECTION  II 


GENERAL  INTEGRAL -EQUAL ION  FORMULATION 


inc 

Let  there  be  a static  magnetic  field  II  ‘ (r)  in  free  space.  Suppose 
a finite  solid  of  uniform  permeability  u is  introduced  into  this  field.  We 
wish  to  calculate  the  induced  field  H^'^^(r). 

l.et  the  space  outside  the  permeable  solid  be  designated  as  region  1 
and  that  inside  as  region  2.  Then  on  the  suriace  S of  the  solid,  the  total 
magnetic  field 


(1) 


satisfies  the  two  boundary  conditions  for  the  tangential  and  normal  com- 


ponents: 


It  - 2t  - 


(2) 


n - In  -J  In-  2n  - 


13) 


ind. 


We  introduce  the  induced  magnetostatic  potential  6 ‘ (r)  such  t'.iat 


,,ind  . _ ind,  . 

H ir)  = -Vo  (r) 


,ind. 


The  conditions  (2)  and  (3)  then  imply  that  6 (r)  is  continuous  across  th.e 

boundary  3 but  that  its  normal  derivative  has  a discontinuity.  It  is  well 
known  that  these  behaviors  are  iniierent  in  the  integral  representatn.-n 


Lad  , I C ’’mV' 

’ ‘ 1 WWT 


( ) 


J 


for  r lying  both  inside  and  outside  the  solid.  That  is  to  say,  the  in- 
duced potential  is  considered  to  arise  from  induced  magnetic  charges  of 


surface  densitv  p (r)  on  the  solid's  surface.  Expression  (5)  is  clearly 
m — 

a solution  of  the  Laplace  equation. 

The  induced  magnetic  field  is  given  by  (4)  and  (5).  As  the  field 
point  r approaches  the  surface  3 we  have  the  following  two  well-known 
limits  in  potential  theory  [2J  : 


(6) 


where  n denotes  the  outward  normal  at  the  field  point  r on  S . The 
sjunbol  P in  front  of  the  integral  denotes  the  Cauchy  principal  value 
whereby  we  must  exclude  from  the  integral  an  infinitesimally  small  cir- 
cular disk  centered  at  r . It  is  precisely  the  integration  over  this  disk 
that  ffives  rise  to  the  term  ~wP  in  (6).  Substituting  (6)  into  (31  we 
obtain  a Fredholm  integral  equation  of  the  second  kind  for  the  induced 
charge  density  [3]  : 


m 


^n 


p + u 


211 ‘■^‘^(rl 
n — 


' 71 


The  solution  is  closely  related  to  tl'.e  total  normal  masr.etic  'uductior.  on 
3.  Eliminating  the  integral  between  '61  :ind  (71  we  obtain 


-i 

i 

I 
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SECTION  III 


RECTANGULAR  SOLID  IN  A UNIFORM  FIELD 

Let  US  specialize  the  integral  equation  (7)  to  the  case  of  a highly- 
permeable  rectangular  solid  placed  in  a uniform  magnetic  field.  The 
situation  is  depicted  in  Fig.  1.  We  set  up  a rectangular  coordinate  sys- 
tem such  that  the  origin  coincides  with  the  center  of  the  solid  and  that 
the  coordinate  axes  are  parallel  to  its  edges.  Let  the  edges  parallel  to 
the  X,  y and  z directions  be.  respectively,  2a,  2b  and  2c  in  length. 

Hence  the  six  faces  of  the  solid  are  defined  by  x = ra,  y = r b and  z = rc. 
Since  any  uniform  field  can  be  resolved  into  three  independent  orthoeonal 
components,  we  can,  without  loss  of  generality,  take  the  incident  mag- 
netic "field  to  point  in  the  positive  x direction. 

From  the  symmetry  of  the  problem  we  need  only  consider  the  values 
of  on  three  faces  of  the  solid.  It  is  convenient  to  introduce  three 

functions  f,  g and  h such  that 


fix,  y) 

2 = =:  C 

II 

N 

X = r a 

(9) 

h(z,  x) 

II 

H 

cr 

These  functions  are  odd  in  x and  even  in  y and  z.  in  practice,  the 
permeability  u is  usually  made  very  high,  so  th.at  the  factor  (a  - 
(u  + u ) in  (7)  can  in  most  cases  be  reolaced  bv  1.  For  exa.mnle,  if 
a > lOOu^,  this  factor  differs  from  l by  less  than  1 . Furthermore  we 
define  where  denotes  the  unit  -v-ector  in  t’le  x -direction. 

Then  (7)  is  equivalent  to  the  set  of  three  coupled  integral  equations: 


r 

r, 

I 


Note  that  the  principal  value  integral  in  vanishes  when  both  r and  r' 
lie  on  the  same  flat  surface. 

In  practice  we  are  often  interested  only  in  the  total  magnetic  flux 
passing  through  the  midsection  of  the  solid  shown  as  ’■he  shaded  area  in 
Fig.  1.  This  total  flux  is 


(X  = 0) 


1 1 j) 


Since  B is  divergenceless,  each  flux  line  crossing  the  sr.aded  area  must 
eventually  come  out  through  the  faces  of  the  solid  for  x '*  0 . T’.ierefore 
'IB  can  -also  be  written  as 


10 


'2  J dzj  dx  B-(y  .b) 
-c  0 

If  we  introduce  the  average  values 


f 


1 

2ab 


A/ 

0 -b 


b 

dy  f(x,  y) 


1 

4bc 


c 

dz  g(y,  z) 


h = 


I 

2ca 


a 

dx  h<z,  x) 


then  by  (8)  and  (9),  for  lu  » u 

n 


expression  Cl 4)  simplifies  to 


SECTION  IV 


METHOD  OF  AVERAGING  FUNCTIONAL  CORRECTIONS 

We  have  on  hand  a rather  formidable  system  of  coupled  two- 
dimensional  integral  equations  in  (10),  (11)  and  (12).  One  method  to  de- 
rive an  approximate  solution  is  obviously  to  replace  these  equations 
with  a set  of  algebraic  equations  by  applying  a suitable  numerical  inte- 
gration formula.  And  yet,  in  practice,  unless  the  problem  exhibits  a 
large  amount  of  symmetry,  as  in  the  case  of  the  electrostatic  analog  of 
a dielectric  cube  [-1],  we  have  to  set  up  quite  a large  number  of  algebraic 
equations  in  order  to  obtain  reasonable  accuracy.  Their  solution  is 
likely  to  tax  severely  the  capacity  of  all  but  the  largest  com.puters.  This 
undesirable  state  of  affairs  is  especially  true  of  the  inductance  calculations 
in  later  sections,  where  the  incident  magnetic  field  is  not  uniform  but 
rather  produced  by  a current-carrying  thin  wire  wound  around  the  rec- 
tancular  solid.  A purely  numerical  solution  will  be  greatly  disadvantaged 
since  the  field  is  singular  at  the  wire,  and  a high-order  numerical  inte- 
gration formula  must  be  used  to  ensure  good  accuracv. 

In  this  work  we  shall  em.ploy  an  anaLtical  approximation  procedure 
known  as  the  method  of  averaging  functional  corrections  [1],  Essentially 
this  method  consists  of  a self-consistent  iterative  calculation  of  the 
average  values  of  the  unknown  functions  over  t!:e  intervals  of  integrati  'ii. 

-As  such,  it  is  ideally  suited  to  calculating  tite  total  flux  :n  ' 1 r i.  since  -.n 
eacli  order  of  th.e  iteration  the  m.etho'i  yte!  is  dtre.'*!-.  •'•i>  I'.  -’rage  v.tlues 
f,  g and  ii. 


The  method  of  averaging  fnnctional  corrections  is  actually  quite 
simple  in  principle.  In  the  first  approxim.ation  we  replace  f , g and  h 
under  the  integrals  in  (10),  (11)  and  (12)  by  their  average  values  defined 


in  (15).  Specifically  we  let 


, y')  - 

( f 

x'  > 

0 

f(x' 

) 

) — 

(-f 

x'  < 

0 

g(y' 

z')  ^ 

g 

1 h 

x'  > 

0 

h(z ' 

(17) 

l-h 

x’  < 

0 

since  f and  h are  both  odd  in  x'  . Under  these  replacements,  the  three 
equations  (10),  (11)  and  (12)  can  be  subsumed  under  the  matrix  eauation 


F - 


F = F 


(18) 


where 


/ f(x,  y)  ^ 

1 

(T\ 

F = 

g^y.  z) 

1 

F = 

a 

\ h(z,  x)j 

[hi 

( 0 


F = 211;  (19) 

o 1 o 1 
■ 0 ! 


and  !\  is  a 3X3  matrix  with  elements 


a „ 

/ dx'  / d, 


Jc 


0 -t) 


(x-  x')“  - (V-  v')“  - 4c' 


3 2 


I X - x')"  - (V  - v')‘ 


A 

4c“  / 


16 


Substituting  (21)  into  (13),  vve  obtain  the  first  iterated  solution  of  (10),  (11) 
and  (12): 


In  the  present  case  F = F . 

‘ o o 

There  are  systematic  procedures  to  improve  the  solution  by  itera- 
tion. But  in  the  present  work  we  shall  not  go  beyond  the  first  iterated 
solution  (23).  We  can  expect  considerable  local  difference  between  the 
exact  solution  and  (23)  at  various  points  on  the  surface  of  the  rectangular 
solid.  On  the  other  hand  we  are  interested  m.ainly  in  the  total  flux  (16) 
which  is  a global  quantity  and  is  obtained  by  integrating  the  induced  mag- 
netic charge  over  the  solid's  surface.  We  believe  (21)  gives  us  a good 
approximation  to  the  flux.  The  present  situation  may  be  com.pared  with 
that  of  calculating  the  capacitance  of  a conductor  where  a rough  approxi- 
mation to  the  surface  charge  density  often  leads  to  an  accurate  value  for 
the  capacitance.  Our  approximation  scheme  is  akin  to  the  Hartree-Fock 
approximation  in  atomic  structure  calculations  and  the  Vlasov  approxim.a- 
tion  in  plasma  dynamics.  The  idea  is  to  replace  an  'unknown  local  or 
microscopic  quantity  under  an  integral  by  its  global  or  macroscopic 
average.  This  average  is  then  calculated  self- consistently  from,  the  inte- 
gral or  integro-differential  equation  thus  approximated. 

For  a detailed  discussion  on  der-.ving  higher-order  approximations 
as  well  as  on  all  other  aspects  of  the  method  of  averaging  functional 


corrections,  the  reader  is  referred  to  the  monograph  of  Luchka[l].  In 
this  book  the  method  is  attributed  to  Yu.  D.  Sokolov. 
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SECTION  V 


NUMERICAL  RESULTS  OF  FLUX  CALCULATIONS 

The  flux  'I'  can  be  obtained  from  (16)  and  (21)  by  working  out  the 
averaged  matrix  K.  Its  nine  elements  are  really  four-dimensional  inte- 
grals. Let  us  first  evaluate  the  set  of  nine  Lvo-dimensional  integrals  K 
in  (20).  It  is  clear  that  this  is  an  enterprise  that  takes  anybody  no  small 
amount  of  effort.  We  shall  not  write  out  the  results  as  they  are  exceed- 
ingly lengthy,  but  shall  only  note  that  they  can  be  expressed  in  term.s  of 
the  following  four  types  of  integrals: 


/dx.  / 

0 -b 


b 

dv' 


C 


0 0 0 
Cx  - x')~  - (y  - y')“  -r  C“ 


3/: 


= tan 


-1 


(x  - a)i  y - b) 


2 0 0 
Cl'x  - a)“  * ly  - b)“  - C" 


172 


- tan 


- 1 


x(  V - b) 


O 0 0 

C|x“  - (v  - b)“  -r-  C” 


1 


- tan 


(x  - a)(  V - b) 


C|(x  - a)“  -r  ( V - b)“  C“ 


U2 


tan 


- 1 


x(v  - b) 


C x”  - (v  - b)”  - C 


1/2 


r 7 

I dy'  I dz' 

■b  - c 


Z - Z ’ 


A“  * (y  - y')“  {z  - z')' 


.3  2 


= fn 


0 

- c)“ 

1 '2 

- b 

o 

( y - b)“  - 

1/2 

(z  c)~  - y - b 

o 

1/2 

1 

0 

1 

O 1 * “ 

^c)-j 

- y - b 

1'- 

(y-  b)“- 

1-.  c)'|  - y-  b 

p 


>.'ext  we  must  integrate  these  complicated  egressions  over  the  solid’s 
surface  to  calculate  K,  in  accordance  with  (22).  The  integrals  can  proba- 
bly be  done  analytically,  but  the  amount  of  algebra  involved  must  be  truly 
enormous.  Since  the  expressions  (24)  have  only  logarithmic  singularities 
at  the  edges  of  the  solid,  we  decide  to  perform  the  remaining  two- 
dimensional  integrations  numerically. 

For  a case  of  immediate  practical  interest  f we  choose  a = b and 
define  a thickness-to- length  ratio 


K = 


£ 

a 


(25' 
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1 

I 


! 

I 


The  values  of  K for  various  < are  shown  in  Table  1.  Note  that 
negative.  Then  f,  g and  h can  be  obtained  by  solving  the  set  of  three 
linear  algebraic  equations  (21).  The  solutions  are  show.n  in  Table  2 to- 
gether with  the  flux  enhancement  factor  $/$  , where 

o 

$ = 4bca  H (26) 

o o o 


is  the  total  incident  flux  through  the  shaded  area  in  Fig.  1 in  the  absence 
of  the  rectangular  solid.  The  ratio  ’5’ also  plotted  versus  < in 
Figs.  2 and  3.  We  see  that  there  is  a dramatic  flux  enhancement  for  a 
thin  rectangular  slab.  But  the  divergence  of  at  sm.all  ■<  can  only  be 

logarithmic.  This  becom.es  evndent  when  we  plot  out  the  quantity 
in  Fig.  4 and  find  that  it  approaches  0 with  <.  Since,  by  (25)  and  (26), 

K.  = — ^5-2 a = b ( 2 1 

■^o  4a"'y  H 

o o 


2 

it  is  reallv  the  total  flu.x  in  units  of  4a“q  H . 

o o 

For  the  case  of  a cube  {<  =1),  we  obtain  a value  of  3.  57  for  _ 

o 

This  may  be  compared  with  the  e.xact  value  of  3 for  an  infinitely  perme- 
able sphere.  The  value  for  a cube  is  bigger  probably  because  of  the  edges 
and  corners.  We  have  also  applied  the  method  of  averaging  functional  cor- 
rections to  the  integral  equation  for  a sphere.  Comparing  the  first 
iterated  solution  with  the  e.xact  result,  we  find  that  the  total  ilu-x  is  smaller 
than  the  exact  value  by  about  8.  3%.  We  hazard  a guess  that  our  flux  re- 
sults in  Table  2 are  sm.aller  than  the  exact  values  by  about  10'”. 
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SECTIOX  VI 


RECTANGULAR  SOLID  IX  THE  FIELD  OF  A 
RECTANGULAR  CURRENT  LOOP 


Consider  a thin-wire  current  loop  wound  tightly  around  the  midsec- 
tion of  the  rectangular  infinitely-permeable  solid.  In  Fig.  1 it  is  to  be 
represented  by  the  perimeter  of  the  shaded  area.  We  want  to  calculate 
the  self- inductance  of  such  a current  loop.  As  is  well  known,  our  task 
amounts  to  calculating  the  total  magnetic  flu.x  passing  through  the  loop  at 
unit  current. 

For  a rectangular  loop  in  free  space  hang  in  the  y-z  plane,  defined 
by  y = ±b'  and  z = =c',  and  carrying  a current  I,  the  vector  potential 
is  easily  found  to  be 


A = A e * A e 

— V— V z — z 


(28) 


with 


it  I 

Ay(x,  y,  z)  = ^ Cn 


o o o 

lx""  (V  - b')“  -r  iZ  c')~ 


V - b' 


o ■>  o 

I x”  - ( V - b')“  - (z  - c')“ 
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V - b ' 


O ■')  1 

c“  ( V - b')"  - (z  - c')“ 


,1 
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V - b' 
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o O 0 

x”  (v  - b')'  - (z  - cM” 


,1  2 


- b' 


U I 

A (x,  V,  z)  = — — in 
z " 4* 


( V - b ' ) “ - < z - c ' ) ■ 


(v-b')"  - (z  c')' 


,1  12 
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It  is  clear  that  for  an  infinitely  thin  wire,  A diverges  logarithmically  at 
the  wire.  In  order  for  (29)  to  represent  the  vector  potential  of  a loop  of 
small  but  finite  wire  radius  R and  wound  tightly  around  a rectangular 
solid  of  cross-section  b by  c,  we  set 


R=b'-b=c'-c 


The  magnetic  field  due  to  such  a loop  in  free  space  is  the  incident  field  in 


the  present  problem,  and  is  given  by 


■)  = J_  V X .\(r) 

— - u 

o 


It  is  clear  that  this  problem  possesses  the  same  type  of  symmetry 
as  the  case  of  a uniform  incident  field  treated  in  Section  III,  namelv,  the 
induced  magnetic  charge  density  on  the  solid's  surface  is  aaain  odd  in  x 


and  even  in  y and  z.  The  only  difference  is  in  the  value  of  the  incident 
field.  Therefore  in  calculating  the  first  iterated  solution  by  the  method 
of  averaging  functional  corrections,  we  again  arrive  at  '21'  but  with  F 


given  by 


/:Hr'\ 


2h'"= 


\ ! 

\ y / 


The  averages  are  over  the  faces  of  the  solid; 


2.2 
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1 ijinc  \ 


0 -b 


4^  /dy  fdz  n;^"dO(  = a) 

- o - c 


^mc 
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-c  0 


..me,  , > 

H w = b) 
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Using  (.31)  and  Stokes'  theorem,  we  obtain 


-pine  ^ 1 

z 2abiU 


J dy  |-A^(a,  y,  c)  - A (0.  y,  c)j 


H = wr ( r dz  A ( a,  D,  z)  - f dv  A (a,  v,  c) 

X 2bcM  \ «/  z J ‘ V - 

o \ - c - o 


j-ine  ^ 1 — / ( ^ (0,  b,  z)  - A (a,  b. 

y 2ca/U  ^ I z z 


For  A and  A given  bv  (29)  and  b'  and  c'  bv  (30),  the  integrals  can  be 
y z 

worked  out  analytically.  We  have 


dz  .A  (a,  b.  z)  — W(a,  b,  c,  R) 

Z 4 T 


y dy  A^m.  y. 


c)  = - W(  a,  c,  b,  R) 
4t 


where 
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The  self-  inductance  in  free  space  is  the  total  flux  through  the  loop  at 
unit  current  and  in  the  absence  of  the  rectanaular  solid.  We  obtain 


u /*b  /.c 

■-o  = r /'»■/'*-  C 

-b  - c 


(X  - 0) 


HI 


-h 


dz  A (0,  b,  z)  - J dv  A (0,  v,  c) 
z - V " 

-b  / 


L = T^[W(0,  b,  c,K)  - W(0,  c.  b,  R)J 

O 2' 


(2n  the  other  hand,  by  (16),  the  self-inductance  L with  the  permeable  core  i; 
4abu 
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SECTION  VII 


^'U-MERICAL  RESULTS  OF  INDUCTANCE  CALCULATIONS 

For  a case  of  immediate  application  we  take  a = b and  R = a/900, 
and  calculate  the  averaged  vector  F in  (32)  and  L in  (38)  numericallv 
as  a function  of  the  thickness-to- length  ratio  < = c/a.  The  results  are 
shown  in  Table  3.  We  then  solve  for  f,  g and  h from  (21)  and  (32'.  The 
matrix  K is  as  given  in  Table  1.  The  self- inductance  L in  (39)  becom^es 
in  the  present  case 
2 

4a”w  _ _ _ 

L = — j — 2[f  <(g  h)]  a = b (401 

The  numerical  results  are  shown  in  Table  3.  We  also  olot  L in  Fig.  5 

o 

and  L in  Figs.  6 and  7. 

The  curve  for  L exhibits  a broad  minimum  at  < 0.  3 and  a sharp 

rise  as  < — 0 . This  behavior  is  at  first  puzzling  since  we  more  or  less 
expect  on  intuitive  grounds  a m.onotonic  dependence  on  < . However,  it  is 
im.portant  to  realize  that  tlie  effect  of  < on  L arises  from  two  sources- 
one  directly  from  the  gecm.etry  of  the  rectangular  solid,  and  another  indi- 
rectly from  the  incident  field  which  varies  with  tlte  loop  geo.metry.  Tc 
study  the  effect  due  purely  to  the  solid,  we  must  factor  out  that  portion 
due  to  the  incident  field.  .Accordingly  we  calculate  and  plot  the  inductance 

enhancen'.ent  factor  L/L  in  Table  3 ,md  Fig.  3.  We  see  that  this  factor, 

o 

purged  of  the  ‘.'-dependence  from,  the  incident  field,  is  indeed  .a  monotonic 

function  of  <.  It  in  fact  closelv  resembles  the  flux  enhancen-. ant  factor 

?'I>  in  Fig.  2.  This  result  is  re.omrkable  since  for  T’  t'r.e  incident 
o o 
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field  is  purely  in  the  x direction,  whereas  for  it  is,  by  Table  3, 

predominantly  in  the  y and  z directions.  It  means  that  the  flux  amplifi- 
cation property  of  the  rectangular  permeable  solid  depends  mainly  on 
the  total  incident  flux  and  not  so  much  on  the  detailed  configuration  of  the 
incident  field. 

From  Fig.  5 we  see  that  L is  a m.onotonicallv  increasing  function 

o ” ■ 

of  K . If  we  write 


L = 


L ^o 


f41) 


we  can  appreciate  that  L need  not  have  cuiy  monotonic  dependence  on  k 

at  all  since  it  is  the  product  of  a monotonically  decreasing  and  a monotoni- 

cally  increasing  function.  The  sharp  rise  of  L as  k -*  0 is  really  due  to 

a faster  rate  of  increase  of  L/L  than  the  rate  of  decrease  of  L . We 

o o 

must  also  point  out  that  for  < v'ery  close  to  0 , L m.ust  eventually  come 
down  to  0.  This  is  because  in  this  region  drops  sharplv  to  0.  Ir. 
fact  we  have  L ~ < \ fn  k ! while  L /L  ~ U n < ' for  < 0 . This  result 


o o 

is  physically  sound  since  the  self-inductance  of  a current  loop  of  zero 
area  must  be  0 . On  the  other  hand  it  is  meaningless  to  push  the  present 
calculations  down  to  very  small  < . Our  treatment  of  the  current  loop  is 
based  on  the  thin-ware  approximation;  this  breaks  down  for  M ~c,  w’len 
the  thin  wire  must  be  considered  as  a massive  cylinder. 

Finally,  what  we  have  calculated  is  actuallv  th.e  external  self- 


inductance. 


None  of  our  results  includes  th.e  interna!  sclf-inductar.ee  due 


to  magnetic  flux  in  the  interior  of  the  thin  wire.  For  a cylindrical  thin 


wire  of  a non-rnagnetic  conductor  in  free  space,  the  internal  self- 
inductance amounts  to  about  u^/8r  per  unit  wire  length  [5].  It  is, 
ever,  negligible  in  our  present  problem. 


how- 


SECTION  VIII 


CONCLUSIONS 

The  general  problem  of  the  interaction  of  a permeable  solid  with 
magnetic  fields  can  be  exactly  formulated  as  a Fredholm  integral  equa- 
tion of  the  second  kind.  While  a number  of  numerical  techniques  are 
available  for  solving  this  type  of  integral  equation,  the  method  of  averag- 
ing functional  corrections  is  found  to  be  particularly  simple  and  power- 
ful when  it  comes  to  evaluating  the  total  magnetic  rlux  passing  through  the 
solid.  In  this  way,  the  total  flux  through  a rectangular  ferrite  slab  is 
calculated  when  the  incident  field  is  uniform  and  when  it  is  generated  by 
a rectangular  current  loop.  It  is  found  that  a thin  ferrite  slab  is  a good 
flux  concentrator. 
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